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Abstract. The Crab nebula has proved to be the nearest to 
a standard candle in VHE 7— ray astronomy. Results on the 
gamma ray emission from the nebula at various energies have 
come in the last decade mostly from imaging telescopes. The 
aim of the new Pachmarhi Atmospheric Cerenkov Telescope 
(PACT) array has been to use the temporal and spatial distri- 
bution of Cerenkov photons in distinguishing between proton 
and gamma ray showers. The array, with timing information 
from 175 mirrors, is ideally suited for precise estimation of 
the arrival direction. Preliminary results from the recent data 
taken on Crab has shown a 12cr signal with the flux in good 
agreement with those of the other experiments. 



1 Introduction 

Crab nebula has been the most intensively studied object in 
all High Energy Astrophysics. There were many attempts to 
detect VHE gamma ray emission from the Crab in the 70s 
and early to mid 80s with varying success. Its certain de- 
tection in the late '80s by the Whipple Imaging telescope 
[Weekes (1988)] established the field of VHE 7- ray as- 
tronomy on a firm basis. Since then many groups includ- 
ing Whipple have been able to establish good signals from 
the source. The energy range from 100 GeV to 7000 GeV 
has been covered by the Atmospheric Cerenkov telescopes. 
There have also been detection of Crab in this energy range 
by EAS experiments. Further, the VHE 7— rays from Crab 
have been explained as inverse Compton scattering of rela- 
tivistic electrons off soft neighbouring photons. [Ong (1998)] 



technique, the aim of the new Pachmarhi Atmospheric Cerenkov 
Telescope (PACT) array has been to use the temporal and 
spatial distribution of Cerenkov photons in distinguishing be- 
tween proton and gamma ray showers for increase of sensi- 
tivity. The array consists of 25 telescopes deployed in a field 
of 80m x 100 m area [Bhat (2000)]. Each telescope consists 
of 7 parabolic mirrors, each of diameter 0.9 m. Each mirror 
is looked at by a fast PMT behind a 3 degree circular mask. 
To minimize the attenuation and distortions in the cables, the 
array is divided into 4 sectors, with each sector servicing a 
group of 6 nearby telescopes. The individual mirror rates 
were kept at about 5 KHz. The pulses from the 7 PMTs in 
a telescope were added linearly to form a 'Royal Sum' pulse 
for each telescope. These Royal Sum pulses were discrim- 
inated to give a rate of about 30 KHz. A four fold coinci- 
dence of these Royal Sum pulses generated the event trig- 
ger. A single sector event rate was about 3 to 4 Hz for most 
atmospheric conditions. Both the timing and pulse height 
information for each mirror are recorded in the nearby sec- 
tor electronics hut serving each sector. All the royal sum 
pulses were brought into the control room and their arrival 
times recorded along with the other housekeeping informa- 
tion. With all 4 sectors functioning, the overall event rate 
was about 9 Hz. Monte Carlo estimates of the energy thresh- 
old for one sector is 900 GeV. When at least 2 sectors are 
demanded in coincidence, the threshold is slightly higher. A 
total of 60 hours of data was taken on the Crab nebula and 
about 20 hours on several background regions, but with the 
same declination as Crab. The Crab data ranged from runs in 
Nov 1999 to Nov 2000. 



2 PACT Experiment 

A new atmospheric Cerenkov array to study cosmic sources 
of Very High Energy (VHE) 7 rays has been set up in Pach- 
marhi in central India. While most of the atmospheric Cerenkov 
experiments in the recent years are based on the imaging 
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3 Estimation of Arrival Direction 

The arrival direction has been calculated using the TDC in- 
formation from the so called Royal Sum pulses and individ- 
ual mirror outputs. Several different methods have been used 
for the estimation of the arrival direction [Majumdar (2001)]. 
We discuss here the preliminary results from one method 
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Fig. 1. (a) a distributions for Monte Carlo generated events assum- 
ing HWHM of 1 (left) and 3(right)ns. (b) The ratio of number of 
source to background events as function of the differences in ex- 
pected and observed delays in telescopes. 



and the results from other methods will be presented in the 
conference. The data used is from the more recent runs in 
fall of 2000 and covers 31 hours (15 nights) and 16 hours 
(10 nights) on source and background respectively when two 
southern sectors were completely operational. Further it was 
required that the events be common in both sectors to in- 
crease the number of degrees of freedom and get a longer 
lever arm. The usable event rate was 1.3 Hz. A simple 
Monte Carlo program was written to understand the angu- 
lar distribution capabilities of the array. The arrival times at 
all the 12 telescopes (in the southern half of the array) were 
generated (for a point source at different zenith angles) as- 
suming a TDC distribution with different values for HWHM. 
Fig 1(a) shows the result of angular reconstruction of Monte 
Carlo events assuming HWHM of 1 and 3 ns. For a TDC 
distribution of HWHM equal to 1(3) ns, it was found that 
73% of the events were within 0.25(0.6)° from the assumed 
source direction. Thus, At, the absolute value of the dif- 
ference between the expected and the observed arrival time 
of the Cerenkov photons at the telescope is an indication of 
the offset of the direction of the event wrt the source direc- 
tion. Before the actual angular reconstruction process, all 
events were binned as to the number of telescopes with At 
< 1, 2 and 3 ns. For eg; there were 146492(74904) events in 
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Fig. 2. (a) q distributions for events in data source(upper) and back- 
ground (lower) (b) The excess number of events v/s a. 



the data from source(background) with at least 4 telescopes 
showing At < 3 ns. The Ratio R, of Source to background 
events where At, the delay difference in at least 4 telescopes 
was < 3 ns is 1.96±0.1. Fig 1(b) shows this ratio R vs the 
delay difference. It is clearly seen that there is a considerable 
excess of source events when At is < 1 ns. The excess, as- 
suming the normalization from the data at < 3 ns, amounts 
to4244±381events. 

The standard angular reconstruction program was used to 
find a, the space angle difference between the events and the 
source direction. Fig 2(a) shows the a distribution for both 
on Crab and off Crab (normalized to total number of events 
on source) directions. Except for differences at small a, the 
distributions are quite similar. It should be noted that these 
distributions are broader than those of Monte Carlo shown 
in Fig 1(a). While the Monte Carlo calculations needed all 
the 12 telescopes, the data analysis demanded at least 4 tele- 
scopes including at least one at either extreme of the ar- 
ray to have At < 1 ns. Fig 2(b) shows the excess events 
(source - background) as a function of a. The total excess 
within 0.6° amounts to 4428 ± 353 events, a 12cr effect. The 
minutely rate is 2.4 ± 0.2. When trigger from at least 2 sec- 
tors are demanded, the energy threshold increases to 1 100 
GeV. The excess from the Crab direction translates to a flux 
of 3.2 ± 0.3 x 10- 1 1 phsec- 1 cm- 2 at >1 100 GeV. 
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4 Discussions 

The VHE 7- ray flux at >1 TeV is 2.1 ± 0.2 ± 0.3 x 10~ n 
sec^cm -2 according to the fit given by the Whipple exper- 
iment. [Hillas (1998)] Given the preliminary nature of the 
above results from the Pachmarhi array, the derived flux for 
the Crab is quite encouraging. However, one should note that 
the PACT array is capable of better angular resolution and the 
methods being developed should increase the sensitivity and 
get more precise value for the flux. Other methods, for eg, us- 
ing the timing jitter [Chitnis (1999)],[Chitnis (2001)], lateral 
distribution parameters etc, are also being evolved. Single 
sector events would help go down in energy threshold. The 
trigger rate with the whole array ( with all the 4 sectors oper- 
ating) is about 50% more than with just 2 sectors. Thus the 
data to be taken in the coming seasons should be able to see 
a strong signal from Crab in a much shorter time. 
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